Massless Dirac fermions in a laser field as a counterpart of graphene superlattices. 



Sergey E. Savel'ev 1 and Alexandre S. Alexandrov 1 
1 Department of Physics, Loughborough University, Loughborough LE11 3TU, United Kingdom 

We establish an analogy between spectra of Dirac fermions in laser fields and an electron spectrum 
of graphene superlattices formed by static ID periodic potentials. The general relations between a 
laser- controlled spectrum where electron momentum depends on the quasi-energy and a superlattice 
mini-band spectrum in graphene are derived. As an example we consider two spectra generated by 
a pulsed laser and by a step-like electrostatic potential. We also calculate the graphene excitation 
spectrum in continuous strong laser fields in the resonance approximation for linear and circular 
polarizations and show that circular polarized laser fields cannot be reduced to any graphene elec- 
trostatic superlattice. Some physical phenomena related to the peculiar graphene energy spectrum 
in the strong electromagnetic field are discussed. 
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A huge surge of interest to graphene (see e.g., [l], Q) 
as the only two dimensional material known so far stim- 
ulates studies of one dimensional and two dimensional 
graphene superlattices (see e.g., [3]). In analogy to usual 
semiconducting superlattices, it is commonly accepted 
that such graphene superlattices should allow to manip- 
ulate the electron spectra and transport properties in 
graphene-based electronics. However, there is still a lim- 
ited number of studies on how time-dependent electric 
field (for instance laser field) can affect both the electron 
spectrum and transport in graphene. Most studies done 
so far are focused on a perturbative response of graphene 
in weak (probe) electromagnetic time-dependent fields 
A growing experimental demand [5] on the data 
analysis of optical properties and electron transport in 
graphene in strong laser fields is the main motivation of 
this work. 

Another motivation is to study temporal-spatial sym- 
metry for massless Dirac fermions in periodic time or 
spatial electric fields. Spatial periodic potentials U(x) 
with its period much larger than the interatomic distance 
transforms the Dirac cone energy spectrum s(p x ,p y ) with 
electron momentum p = (p x ,Py) into a set of minibands 
£n(Px,Py) [3] where the quasi- momentum p x = ftk x or 
the electron wave vector k x takes values within the first 
superlattice Brillouin zone —tt/L < k x < tt/L, while the 
y-momentum p y = ftk y is unbound. Here L is the spatial 
period of an electric field E x (x) and an integer n numer- 
ates different subbands (which are usually separated by 
gaps). In contrast to the above picture, applying time pe- 
riodic electric fields E x (t) should change the spectrum of 
electrons in accordance with the Floquet theory stating 
that the electron energy should become a quasi-energy e 
bounded within its Brillouin zone —tt/T < £ < tt/T with 
£ = e/h and T is the temporal period of homogeneous 
electric field E x (t) applied along the x-axis. Therefore, 
the electron spectrum in time-periodic laser fields could 
be written in the form k x = k x (s,k y ). This analogy raises 
several questions, including the possibility of gaps in the 
spectrum for either momentum k x or quasi-energy e. 



In this context, we would like to note that a study 
of conventional gapped semiconductors in strong laser 
fields has a long history. It was shown that beyond linear 
response theory, a strong harmonic field opens energy 
gaps caused by the high-frequency Stark effect (e.g., [6(). 
The location of the gaps in momentum space and their 
size can be controlled independently by changing either 
frequency or amplitude of laser fields. 

In this letter we show that two problems of electron 
excitation spectra in graphene in either (i) spatial or (ii) 
temporal periodic potential have a deep analogy allow- 
ing to obtain an unknown spectrum if the dual spec- 
trum is calculated or measured. This opens a pos- 
sibility to indirectly measure the spectrum of compli- 
cated graphene superlattices with no expensive and time- 
consuming nanofabrication techniques needed to gener- 
ate a desirable spatial variation of electric fields on a 
nano-scale but rather applying easily adjustable tempo- 
rally periodic laser electric fields. Then, we consider 
two specific examples: (1) an analytically solvable model 
for periodic piece-wise constant (square wave-like) po- 
tentials, which can be realized by using a sequence of 
laser pulses or by making spatial-potential steps ±E/o, (2) 
monochromatic linearly polarized laser fields (E x (t) oc 
sincjt), which can be analyzed by using the resonance 
approximation (e.g., [6|). As an example of a laser field 
problem irreducible to any spatially periodic electrostatic 
graphene superlattices we consider circularly polarized 
electric fields generating a quasi-isotropic gap in the en- 
ergy spectrum in contrast to an anisotropic gap with gap- 
less points for linearly polarized laser fields. 

Spatial-temporal duality. — Let us consider two related 
problems of graphene subject to either spatially periodic 
E x {x) = —dU/dx or time periodic E x (t) = —dA x /dt elec- 
tric fields. Here we introduce the static electric potential 
U(x) and time-dependent vector potential A x (t) which 
are periodic functions of their arguments, i.e., U (x+L) — 
U(x) and A x (t + T) = A x {t). In the low energy limit, 
the behaviour of the charge carriers can be described by 
the 2D Dirac equation: v F (a(p- A))V + = ihdV/dt 
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with a two component wave function \I> = (^a^b) cor- 
responding to electrons in two triangular graphene sub- 
lattices and vp is the Fermi velocity. For the case of 
spatially periodic electric fields (graphene superlattices), 
i.e., A x = 0, the wave function can be written as fol- 
lows S&a,b = (PA,b(x) exp(— i£t + ikyy) and the coupled 
equations for (fA,B have the form: 



dx 



VA 



ih 



-£ 



y 

U(x) 



-£ + U{x) 



¥B 



(1) 

Here and below we use % = vp = 1 and measure wave 
vectors in the unit reciprocal lattice vector of the super- 
lattice, tt/L, or of the graphene lattice, ir/a. For the 
case of time dependent electric fields E(t) = —dA x /dt 
and U(x) = 0, the solution can be written as ^a,b = 
4>a,b ex.-p(ik x x + ik y y), where 
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Both sets (pQ) and (j2j) of differential equations belong 
to the same class of ordinary differential equations with 
periodic coefficients so that one can apply the Floquet 
(or Bloch) theory to classify all possible states. To 
see a deeper analogy we can use new variables (j)A — 
(<j>A — i(\>B)/\f2 and — {<t>B — i(\>A)/\[2 for time- 
dependent potentials, so that both sets (jTJ [2]) can be 
written in the same form 
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where + C) = /i(£) is a periodic function of £ with 
the period C. Here, £ = x. ; /.i = U, Q = —£, f3 = k y , and 
C = L for spatial-dependent electric fields while £ = £, 
\i = — A x , Q = k x , j3 = —ik yi and £ = T for time- 
dependent homogeneous electric field. All the solutions 
of equations ((3]) can be classified using Floquet exponents 
+ = wp(iK C) f a,b{£>) with Floquet multipliers 
—tt/C < K < tt/C. Spectrum linking Floquet multiplier 
(quasi-momentum k x for spatially varying potentials and 
quasi-energy £ for time-dependent laser fields) depends 
on the shape of the periodic potential \i, but, in gen- 
eral, can be written as F(K,Q,f3) = 0. Therefore, the 
spectrum of Dirac electrons in the periodic superlattices 
with U = P(x) and the spectrum of Dirac electrons in 
laser fields with A x = —P(t) described by the same func- 
tion P are defined by the same spectral equation, i.e., 
F{k x ,—£,k y ) = for spatial oscillating potentials and 
F(—£,k x ,ik y ) = for laser fields. For instance, if we 
measure a quasi-energy spectrum for a laser field having 
a certain time-profile P(t), then we can easily obtain the 
spectrum of graphene superlattices with a spatial poten- 
tial of the same shape U(x) = —P{x). 

Periodic spatial step-like potential versus pulsating 
laser fields. — Generally it is not possible to write down 
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FIG. 1: (Color online) (a) Mini-band electron structure for 
a graphene sample in the step- like spatial potential U(x) — 
Uo sign (sin (2-kx/L)) with Uo = 1, L = 1, k y = 1. The 
gaps (Si, S2, S3, ...) between different zones are clearly seen, 
(b) Reduced zone momentum- quasi- energy spectrum k x (£) 
for graphene in the laser field A — Aosign(sin(27r£/T)) with 
Ao — 1,T — l,k y — 1. There are no gaps for momentum 
but not all values of quasi-energy can be reached resulting in 
gaps in quasi-energy Ai, A3, A3 marked on the figure. The 
physical meaning of the energy gaps are readily seen in the 
extended zone representation (c). 



the solution of Eqs. §3§ for k y 7^ and, thus, to de- 
rive the spectral equation F(K,Q,j3) = for an arbi- 
trary function However, there are several partic- 
ular cases when we can obtain analytical solutions and, 
thus, the corresponding spectra of Eqs. (|3]). One in- 
teresting example is the step-like spatial periodic func- 
tion U(x) = /7osign[sin(27rx/L)] and the dual periodic 
temporal potential A x = Aq sign [sin (27rt/T)], both re- 
sulting in = /iosign[sin(27r£/£)] (with /xq = Uo 
or fio = Ao). Here sign(x) = 1 for x > and 
— 1 for x < 0. A solution for /a,b can be written 

as f A ,B 
r (±) 



= ^A^B;! ex P(^^lO + C A,B-A ex P(-^l)£ witn 

(T*i - iP)c£l/(Q - fio) for -C/2 < £ < 



and /a,B;2 
,(±) 



Ca^b-,2 exp(m20 + ^a,B;2 exp(-m 2 )£ with 
= (=f^i - if3)C^/(Q + fjL ) for < £ < C/2. Here 

we introduce K\ 2 as K\ 2 = V(Q ± Vo) 2 ~ P 2 where the 
sign "— " should be taken for subindex 1 and for 
subindex 2, while subindexes 1 and 2 refer to the regions 
-C/2 <£<0and0<£< C/2 respectively. It is im- 
portant to stress here that can become imaginary 
for spatially oscillating fields when f3 2 = k 2 > (Q±Eo) 2 , 
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but ^i 5 2 are always real for time-oscillating fields since 
f3 2 — (iky) 2 = —k 2 . Using both, the matching conditions 
for /a,b ; i,2 at £ = and the quasiperiodic conditions 
f AtB (i = C/2) = exp{iKC)f A ,B(£ = -C/2) we can de- 
rive four homogeneous equations for four variables C^. x 2 
which can be solved only if its determinant is equal to 
zero, which defines the spectral equation: 
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Following our duality rule we have to substitute K = 
kxiQ — —£->fi — k y and C = L for graphene super lat- 
tices (-E(x)-electric fields) and K = — £, Q = k x , j3 = ik y 
and C = T for time-dependent laser pulse fields. As 
mentioned above, can become imaginary for spatial 
oscillating fields, which results in the right-hand-side of 
equation (|4]) to exceed 1, thus, giving no roots to this 
equation for such values of energy (energy gaps Si , 62 , . . , 
see Fig. la). In contrast, the right-hand side of equation 
(|4|) never exceeds 1 (moreover is always less than 1 for 
j3 2 = — k 2 < 0), thus, the equation has solutions for any 
Q = k x and k y producing no gaps in the momentum. 
However, some of the values of quasi-energy £ = — K 
near the boundaries ±T/7r of the Brillouin zone cannot 
be reached resulting in gaps in quasi-energy (see Fig. lb). 
Note that the energy or quasi-energy spectrum of equa- 
tion (j4]) has no gaps for either spatial or temporal oscil- 
lating electric fields if k y = 0. This last property could 
be obtained for any spatial U (x) and temporal A x (t) de- 
pendence of electric fields, which can be proven by using 
general solution (for k y = 0) obtained in [7[ for the Dirac 
equitation (j2j). Sometimes it is more convenient to re- 
place the quasi-energy spectrum (Fig. 1) in the reduced 
zone representation —T/ir < e < T/tt by the spectrum in 
the extended zone representation shown in Fig. lc where 
energy gaps Ai, A2, ... are clearly seen. 

Harmonic linearly polarized laser field. — Unfortu- 
nately, an exact solution for the harmonic poten- 
tial A x (t) = Aocos(ut) with frequency uo = 2n/T 
or its dual spatially periodic potential U(x) = 
Uo cos(2ttx/L) is unknown. However, we can use the 
so-called resonance approximation to calculate the spec- 
trum of Eqs. (j2j). Namely, we can seek a solu- 
tion in the form (f)A,B = exp(i£t)[D^ B exp(iuot / 2) + 

D { ~1 exp(-iut/2)}. When ignoring all fast oscillating; 
terms [i.e., exp(±3zcj£/2), exp(±5zco>£/2), ...], we can ob- 
tain a closed set of four homogeneous equations for four 
variables D^bi which can be solved only if its determi- 
nant is equal to zero resulting in the spectrum: 
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where ± in front and inside the brackets should be taken 
independently, producing four independent branches 
shown in Fig. 2b, c. 

In agreement with our general statement, the oscillat- 
ing laser field opens a gap in the spectrum, which can 
be controlled by both the frequency uo and the field in- 
tensity Aq. The gap between different branches of the 
spectrum drops to zero for k y = and any field inten- 
sity Ao(as expected from our general consideration) (Fig. 
2b), while gaps are seen for nonzero k y (Fig. 2c). Phys- 
ical meaning of four different branches can be easily un- 
derstood in the limit of Aq — >> 0. Indeed, in addition to 
the usual energy cone (Fig. 2a(left); also, the dark-blue 
line for electrons and the dark-green line for holes in Fig. 
2a(right)), we have to consider the energy spectrum of a 
hole and a photon (a cone shifted up by uj shown in Fig 
2a(right) in light-green) and the spectrum of an electron 
and an emitted photon (a cone shifted down by uo shown 
in Fig. 2a(right) in light-blue). Gaps can be formed 
at the crossing points of these spectra for nonzero field 
intensity (so-called the high-frequency Stark effect). In 
this case, we recover the spectrum shown on Fig. 2b for 
the energies £±uo/2 since the wave function contains the 
terms D A B (±) exp(i(£ ± uo/2)t). 

Circular laser fields. — An interesting example of a 
problem, when a time-dependent laser field has no 
analogue in a spatially periodic electrostatic field, is 
the circular polarized harmonic laser field described by 
a vector potential A = (A x (t), A y (t)) with A x (t) = 
(A0/V2) cos uot and A y (t) = (A /V2) cos(uot + a). Seek- 
ing a wave function in the form ^a,b = 4>A,b exp(ik x x + 
ikyy) we recover equation (j2]) with replacement k y by 
k y — A y (t). Using ())a,b = exp(i£t)[D^ B exp(iuot/2) + 

D^a b ex P( — ioot/2)] and ignoring higher harmonics in the 
resonance approximation, we again obtain the set of four 
homogeneous equations for D^ B with the spectrum 
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Note that the last spectrum has no more requirements of 
zero gaps for k y = 0, so the four quasi-energy branches 
can be well separated for a properly chosen phase shift 
a (see Fig. 3). Thus, the spectrum becomes even more 
tuneable by changing frequency, amplitude and/or phase 
of the circular polarized laser field. 

Perspectives and conclusions. — The optical gap in the 
spectrum of graphene electrons generated by the elec- 
tromagnetic field can be responsible for unusual electron 
dynamics and kinetics, for instance, the pseudo-periodic 
£(k x (i)) Bloch oscillations can be induced by applying a 
small DC electric field. Stationary distributions of elec- 
trons in strong laser fields can be very peculiar compared 
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FIG. 2: (Color online) (a, left) The sketch of the excitation 
of the electron Dirac cone in a harmonic laser field with the 
frequency cj, arrow shows electron transitions under the elec- 
tric field, (a, right) Crossing of electron-photon terms under 
the action of the electric field: the Dirac cone shifts both 
up (representing the hole and photon state) and down (rep- 
resenting an electron which has already emitted a photon), 
the crossing of all these three cones [the original Dirac cone 
(in dark-green and dark-blue) and the two cones shifted by uj 
(in light-green and light-blue)] results in the opening of gaps 
between the corresponding energy bands, (b) Electron spec- 
trum 6{k x ) for graphene in a linearly polarized electric laser 
field with o; = l,Ao = l,fc y =0 (all measured in tt/T): there 
are no gaps between the energy zones in agreement with our 
spatial-temporal symmetry conclusions, (b) The same as (a) 
but for nonzero y- momentum k y = 1, resulting in the opening 
of gaps between the energy bands, (c) The same as (a) with 
zero y-momentum, but for the circular polarized laser field 
a — 7r/2 also resulting in the opening of the gaps between 
zones since there is no analogy between the circular polarized 
light and any periodic electrostatic potential. 

to the electron distribution [6] in conventional (gapped) 
semiconductors. In the latter case with a large band- 
gap one can readily reach a quasi-equilibrium Fermi- 
Dirac distribution of electrons pumped into the conduc- 
tion band and holes in the valence band, and an insulat- 
ing state due to the isotropic optical gap (the so-called 
"optical" insulator [8[). However, such degenerate distri- 
butions seem to be improbable in graphene, since pho- 
toelectrons can recombine with photoholes with about 
the same rate as their energy relaxation rate. More- 
over, the spectrum predicted here can be directly mea- 
sured by applying a weak probe field in additional to 
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strong time-periodic laser or spatially periodic electro- 




FIG. 3: (Color online) Electron spectra E(k x ,k y ) for 
graphene in linear polarized (a[upper and bottom energy 
branches] ,b [two middle energy branches]) and circular polar- 
ized (c[upper and bottom energy branches] ,d [two middle en- 
ergy branches]) laser fields with parameters Ao = 0.03, uj — 
0.1 (measured in units hvFTv/a with a being the atomic lattice 
constant). Two points for k y = where two middle energy 
branches touch one another are clearly seen for linearly po- 
larized electromagnetic fields, while these branches are well 
separated for any k x , k y for the circular polarized electromag- 
netic fields (a = 7r/2) (c,d). 



static fields. The duality between the electron spectrum 
in spatially periodic electrostatic fields and time-periodic 
laser fields provides an experimental access to the com- 
plex spatially periodic graphene superlattices that was 
impossible before due to the expensive, hardly controlled 
and not easily tuneable nano- fabrication procedure. 
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